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Abstract. We present a new EROS-2 measurement of the microlensing optical depth toward the Galactic Bulge. Light curves 
of 5.6 x 10 6 clump-giant stars distributed over 66deg 2 of the Bulge were monitored during seven Bulge seasons. 120 events 
were found with apparent amplifications greater than 1.6 and Einstein radius crossing times in the range 5d < t E < 400 d. This 

■ is the largest existing sample of clump-giant events and the first to include northern Galactic fields. In the Galactic latitude 
range 1.4° < \b\ < 7.0°, we find t/10~ 6 = (1.62 ± 0.23) exp[-a(|£>| - 3deg)] with a = (0.43 ± 0.16)deg-'. These results are 
in good agreement with our previous measurement, with recent measurements of the MACHO and OGLE-II groups, and with 

'_h ' predictions of Bulge models. 

■ Key words. Galaxy:bar - Galaxy: stellar contents - Galaxy:structure - Cosmology:gravitational lensing 
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' 1 . Introduction events, the lensing object passes near the line of sight towards 
• the background star, causing a transient magnification of the 
Gravitational microlensing of stars is an important tool for star ' s p r i ma ry image as well as creating a secondary image, 
constraining the quantity and characteristics of faint compact At Galactic scales, neither the image separation nor the im- 
objects between the stars and the observer. In microlensing age size are norma ii y resolvable, so the only easily observable 
effect during a microlensing event is an apparent transient am- 

* Based on observations made with the MARLY telescope at the plification of the star's flux. The characteristic timescale giving 
European Southern Observatory, La Silla, Chile. the effective duration of the amplification is proportional to the 
** Now at Instituut voor Sterrenkunde, Celestijnenlaan 200 B, B- square root of the lens mass and inversely proportional to the 
3001 Leuven, Belgium transverse velocity. For Galactic stellar lenses, the timescales 
*** Now at Research School of Astronomy and Astrophysics, on me orc [ er G f a few weeks 

Australian National University, Mount Stromlo Obs., Cotter Rd„ Microlensing searches were originally proposed 

Weston, ACT 261 1, Australia JPaczvnskill98fj|) as a tool for detecting dark matter in galactic 

I Now at Max-Planck-Institut fur Astronomie, Koenigstuhl 17,D- _ , .,»„,,- 
.„,,,., ., „ _ halos. Searches for the lensing of stars in the Magellanic 
691 17 Heidelberg, Germany I a ii 1 ° 

* Now at Max Planck-Inst.tut fur Astronomie, Koenigstuhl 17,D- Clouds by the MACHO <Alcock et aljEig) and EROS-2 
69117 Heidelberg, Germany dAfonso et alJ 12003d: iTisserand et alJ 1200 6) projects have 

* Now at Division of Medical Imaging Physics, Johns Hopkins P laced constraints on the fraction of the Milky Way halo that 
University Baltimore, MD 21287-0859, USA can be comprised of faint compact objects. Searches for objects 

II Now at San Pedro de Atacama Celestial Exploration, Casilla 21, in the halo of M31 are also being perf ormed with candidat e 
San Pedro de Atacama, Chile events reported by the VATT-Columbia ( Ugle sich et alJl2004l) . 
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WeCA PP dRiffeseretalJl2003l). AG APE dCalchi Novafi et al 



2005I) . MEGA dde Jong et al J 12006) and Nainital dJoshi et al 



2005) collaborations. The AGAPE and MEGA collaborations 



presented efficiency calculations allowing them to constrain 
the content of the M31 and Milky Way halos. 

Microlensing has also been used as a tool to investigate 
compact objects in the visible regions of the Milky Way. 
Besides the Galactic Bulge that is the subject of this paper, mi- 
crolensing of stars in the spira l arms of the Galaxy has also 
been studied ( Deru e et al. 

Microl ensing surveys toward the Gala ctic Bulge were first 
proposed dPaczvnskilll99ll lOriestl 1 1 99lh as a probe of ordi- 
nary stars in the Galactic Disk, though it was soon realized 
dKiraga & Paczvriskill 19941) that lensing by stars in the Bulge 
itself is of comparable importance. The optical depth, i.e. the 
probability that at a given time a star at a distance D s is ampli- 
fied by more than a factor 1.34 is 
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p(x)x( 1 - x) dx 
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where p is the mass density of lenses and x = D\/D s is the ratio 
between the lens and source distances. Very qualitatively, the 
integral for the Disk contribution to the optical depth is 

GM disk 



Tdisk 



(2) 



C 2 /ldisk 

where Md; s k and /idisk are the total mass and scale height of 
the Disk lensing population. Lenses in the Bulge are near the 
source star so the Bulge contribution is 
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where Mb u i ge , Pbuige and R\, u i g e are the mass, density and line 
of sight thickness of the Bulge lensing population. For Bulge 

stars near the direction of the Baad e window (i = \°,b = 

J r 1 

3.9°), the most recent calculations (Evans & Belokurov 2002; 



Biss antz & GerharJ2002HHan & Goul d 2003; Wood and Mao 
2005I) give total optical depths in the range 1 < t/10~ 6 < 2, 
with about 60% of the rate being due to lensing by Bulge stars 
and the remainder by Disk stars. The expected optical depth 
naturally has a strong dependence on the Galactic latitude, b, 
falling typically from r ~ 5 x 10~ 6 at (I ~ 0, \b\ = 1°) to 
t ~ 5 x 10~ 7 at (I ~ 0, \b\ = 6°). 

The first measurements of the optical dept h towards the 
Galact i c Bulge by the OGLE co llaboration dUdalski et aid 
Il994bi l2000t IWozniak etlil l200ll) . the MACHO collabora- 
tion dAlcocketal]ll997ll2000al) and the MOA collaboration 
dSumi et alTT2003l) yielded optical depths significantly higher 
than these estimates, suggesting that a fundamental r evisio n 
of Galactic models might be necessary dBinnev et alJ feoOO). 
However, the interpretation of the first results was difficult be- 
cause of two effects. First, stars in the direction of the Galactic 
Bulge are not necessarily located in the Bulge itself so fore- 
ground and background contamination must be taken into ac- 
count. Second, and more importantly, photometry of the nu- 
merous faint stars in the crowded Galactic Bulge fields is com- 
plicated by "blending" where the reconstructed stellar flux re- 
ceives contributions from more than one star. This makes the 



effective number of monitored stars greater than the number 
of cataloged stars and generates a complicated relationship be- 
tween real and reconstructed amplifications. 

These problems can be largely avoided by considering only 
microlensing of clump giant stars. Such stars are identified 
by their well-defined position in the color-magnitude diagram. 
This position ensures that they are most likely in the Galactic 
Bulge. Their large flux also makes blending problems relatively 
unimportant. 

The first published optical depth using clump giants was 
based on only 13 events observed by the MACHO group 
dAlcock etalJ ll997h . The optical depth, 3.9+J* x 10~ 6 at 
(l,b) = (2.55°, -3.64°), was still considerably higher than 
expectations. Since then, measurements using clump giants 
have tended to give results in better agreement with models- 
Based on 16 events, the EROS-2 collaboration dAfonso et alJ 
l2003bl) gave an optical depth of 0.9 4 + 0.29 x 10" 6 at (lb ) = 
(2.5°, -4.0°). The MACHO group dPopowski et al.ll2005l) re- 
cently presented their results using 62 events with clump gi- 
ant sources. Their optical depth is t = 2.17^'^ x 10~ 6 at 
(I, b) = (1.5°, -2.68°), in good agreeme nt with ca lculation s an d 
with their earlier preliminary results dPopowskietalJ EoOO). 
The OGLE-II group dSumi et alJ2006l) recently used 33 events 
to obtain t = 2.55+°;^ X 10~ 6 at (l,b) = (1.16°, -2.75°), again 
in agreement with models. 

In this paper we present a measurement of the optical depth 
for microlensing toward the Galactic Bulge based on the to- 
tality of the EROS-2 clu mp-giant da t a. Mo re details on this 
analysis can be found in Hamadache (2004). The 120 events 
used comprise the largest sample of clump-giant events studied 
to date. Section 13 describes the data collection and process- 
ing required to produce the light curves. Section [3] describes 
the definition of the clump-giant sample to be used to mea- 
sure the optical depth. Section [4] lists the criteria used to iso- 
late a sample of microlensing candidates. Section [5] describes 
the calculation of the detection efficiency, the optical depth and 
the timescale distribution. Section [6] discusses the possibility 
that stellar blending may affect the results. Finally, Section Q 
compares our results with those of other groups and with the 
predictions of Galactic models. 



2. Data 

The data were acquired by EROS-2 with the 1 m MARLY tele- 
scope at La Silla, Chile. The imaging was done si multaneously 
by tw o cameras, using a dichroic beam-splitter dBauer et alJ 
1997). Each camera was composed of a mosaic of eight 
2048x2048 pixel LORAL CCDs, with a pixel size of 0."6 
yielding a field of 0.7°(a) X 1.4°(<5). Each camera observed 
through one of two non-standard, broadband filters, R sms and 
B elos . The filter R el0S covers the range 620-920 nm and has a 
mean wavelength near that of Cousins /, while B eros covers 
the range 420-720 nm and has a mean wavelength between 
Johnson/Cousins V and R. The photometric calibration was ob- 
tained directly for 20% of our fields by matching our star cata- 
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galactic longitude 

Fig. 1. Map of the EROS-2 Bulge fields in Galactic coordinates. 
A total of 66 fields were monitored. The first number in each 
field is the field number and the second is the number of clump 
giants in the field divided by 10 4 . The field 121 with 4 x 10 4 
clump giants is outside the plot at (I, b) = (11.5, 2.5). 



logue s with those of the OGLE-II collaboration ( Udalski ^t alJ 
2002). To a precision of ~ 0.1 mag, we find 



^eros — ^ogle 



B, 



= Vogle-0.4(V-/) ogle 



(4) 



The calibration is sufficiently uniform that it can be extended 
with confidence to the remaining fields. We note, however, that 
none of our conclusions depend on the calibration precision. 

The observations reported here concern 66 fields, moni- 
tored between July 1996 and October 2002. Figure^shows the 
location of the 66 fields in Galactic coordinates, (/, b). Twenty- 
two of the fields are in the north (b > 0) and 44 in the south. The 
corresponding data set contains 6 x 10 7 pairs of light curves, 
of which 5.6 x 10 6 are Bulge clump giants that we use in the 
present analysis. 

The fields are observable at La Silla from mid-February to 
mid-October. An average of one point every four nights was 
taken for each field during the observing season though the 
sampling frequency occasionally reached more than one point 
per night at mid-season when the Galactic Bulge is observable 
throughout the night. 

The image photometry was performed with software spe- 
cifically designed for crowded fields, PEIDA (Photometrie et 
Etude dTmages Destinees a l'Astrophysique) (Ansari 1996). 
The first step consisted of creating reference images for each 
of the 2112 subfields of size ~ 0.035 deg 2 , i.e. one quarter of a 
CCD. A reference image was constructed by adding 15 high- 
quality images of the subfield. Star catalogs for each subfield 
were established using the reference image. Photometry on in- 
dividual images was then performed by imposing the stellar 
positions found on the reference image. 

A first estimate of the uncertainty of individual flux mea- 
surements is that due to the photometric fit taking into ac- 
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Fig. 2. The dispersion, oV/base, of clump-giant light curves 
about the base flux as a function of R aoti and of B e ros- Each 
point represents the median dispersion of clump giants for one 
subfield consisting of a CCD-quadrant. 



count only photon counting statistics. This clearly underesti- 
mates the uncertainty since it ignores other contributions im- 
portant in crowded field photometry, e.g. imperfect knowledge 
of the point-spread-function (PSF). The factor by which the 
uncertainties must be increased is determined by comparing in- 
dividual flux measurements with the flux measurements of the 
same star on the reference image. The fit errors are increased by 
a magnitude and image dependent (but star independent) factor, 
chosen so that for each image the distribution of the difference 
between flux and reference flux is consistent with that expected 
from the renormalized flux uncertainties. By construction, this 
procedure yields, for most stars, a x 2 P er degree of freedom 
near unity for a fit to a time-independent flux. 

The uncertainties on flux measurements are reflected in the 
dispersions of light curves about the baseline flux (estimated 
from the mode of the flux). Figure |2 shows, for each subfield, 
the median dispersion for clump giants in the two colors. The 
photometric precision is typically 2% but degrades consider- 
ably for B eros in highly absorbed fields. We will be concerned 
only with microlensing apparent amplification factors > 1.6 
so the precision in most fields is sufficiently good to yield 
high quality light curves in both colors. The 20% of the mi- 
crolensing candidates with the lowest quality light curves were 
also photometere d with our differential photometry package 
tLe Gu illou 2004J) so as to confirm their microlensing origin. 
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Fig. 3. A typical color-magnitude diagram of a subfield (field 
610, CCD 4 quadrant 2). The cross shows the fitted center of 
the clump and the thin lines show the cuts defining the clump 
giants used in this analysis. 

Examples of light curves containing microlensing events 
are shown in Figs. 1211231 

3. Clump giant selection 

The sample of clump-giants to be studied was selected us- 
ing the color-magnitude diagrams for each CCD-quadrant sub- 
field. Figure[3]shows a diagram of a typical subfield. The cen- 
ter of the clump was determined by modeling the density of 
stars in color-magnitude space as a power law representing a 
smooth background plus a Gaussian representing the clump. 
The Gaussian had the form 

, . ((B - R) ems - a) 2 

density oc exp 



2<rl 



x exp 



(R eIm + 1.9(B-R) ems -B) 2 



Icrl 



where a (cr a ) and B (erg) are the desired clump centers (widths) 
in the space defined by (B - R) ems and the reddening-free mag- 
nitude 7? eros - 1 .9(B - R) eros . The dereddened magnitude is used 
so as to be relatively insensitive to variable reddening across 
the subfields. 

Figure shows the fitted R eros of the clump center as a 
function of the fitted {B - R) aos of the clump center. The dis- 
tribution follows the expected reddening law indicating that 
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Fig. 4. The top panel shows R ems of the clump center vs 
fieros - ^eros of the clump for the 2109 subfields (CCD quad- 
rants). (Three quadrants were not reconstructed.) The distri- 
bution follows the expected reddening law. Quadrants marked 
with a x were eliminated from further consideration. The bot- 
tom panel shows the dereddened magnitude as a function of 
Galactic longitude. The increasing flux with longitude is be- 
lieved to be due to the bar structure since stars with I < are 
farther than stars with I > 0. 



variations of clump position are due to differential reddening. 
The 46 subfields that are outliers in this plot are excluded from 
further consideration leaving 2063 subfields. All 2063 color- 
magnitude diagrams were visually inspected to ensure that rea- 
sonable values for the clump position were found. We note that, 
in spite of the large variation of clump position, the clump is vi- 
sually distinct in each of the 2063 diagrams. The variations of 
clump position lead to varying photometric precisions (Figure 
13 but this is taken into account in the microlensing detection 
efficiency (Section[5}. 

Figure |4j) shows the dereddened clump magnitude as a 
function of Galactic longitude, t. The longitude dependence 
is believed to be caused by the bar structure of the Bulge 
where stars with I > are b righter than stars with / < be- 
cause they are nearer (Stan ek et al.l ll994). The 0.38 mag dif- 
ference over 12° in longitude corresponds to a bar orientation 
of 49° ± 8° with respect to the line of sight to the Galactic 
Center. This angle is in agreement with the original OGLE- 
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Fig. 5. The mean clump color (top panel) and the number of 
clump giants per CCD-quadrant (bottom panel) as a function 
of Galactic latitude for the longitude range -3° < I < 3°. 
The number of giants is symmetric about b — while north- 
ern fields are more reddened than southern fields. The large 
dot marks the position of the subfield containing the Baade 
Window. 



I results lStaneketal.ll 1994 and with the recent infrared star 
counts from GLIMPSE lBeniamin et all2005l) . 

Figure [5^ shows the clump color as a function of Galactic 
latitude. The northern fields are more absorbed and reddened 
than the southern fields. 

The cuts used to define our sample of clump giants, illus- 
trated in Figure |3] are 



4. Event Selection 

After the production of the light curves, microlensing candi- 
dates among the 5.6 x 10 6 clump giants were found by the pro- 
cedures described in this section. After elimination of images 
with instrumental problems or poor seeing, the first step was to 
calculate the mode of the red and blue fluxes that was used as 
a first estimate of the baseline fluxes. 

The light curves were then subjected to a "filter" that selects 
light curves having one or more groups of points with fluxes 
sufficiently above the baseline flux, /b aS e- Such groups were ini- 
tiated with any point i having a flux fi greater than /b ase + 2.5cr, 
where cr, is the estimated uncertainty of /;. The groups ended 
with the first subsequent sub-group of three points all with 
Si < /base + 2cr ; . The relative significance of the groups on a 
light curve was defined by 



LP N = AHog2 - J] log 



erfc 



fi /base 



where N is the number of points in the group. For each red 
or blue light curve, groups having at least 5 points with fi > 
/base + 2.5cr, were then ordered by decreasing significance. A 
light curve was then retained for further analysis if one of the 
two most significant red groups had a temporal overlap of at 
least 20% with one of the two most significant blue groups. 

Roughly 3% of the light curves pass this filter. The filter's 
effective threshold is sufficiently low that, for clump stars, the 
filter passes essentially all microlensing events with apparent 
amplifications greater than 30% occurring during well-sampled 
parts of the observing period. Significant losses of efficiency 
occur due to bad weather, equipment failure and limited visi- 
bility of the Galactic Bulge at the beginning and end of each 
season. 

Light curves passing the filter were then fit in each color 
with the microlensing light curve corresponding to uniform 
motion of an unblended point source and single lens. The total 
flux as a function of time is 



fit) = /ba 



M 2 +2 



base 



u ^u 2 + 4 



where the impact parameter u{t) is 



U 2 (t) = Wq + 



(t - to) 2 



(8) 



(9) 



\R - ^clumpl < 1 . 

\(B -R)-(B- /?) iumpl < 0.25 , 

^-^clump > ~4.[(B-R) - (B -/?) c lump] ~ 1 



(5) 



(6) 



(7) 



where /? c i ump and (B - R) c i ump define the clump center. The 
number of clump giants as a function of Galactic latitude is 
shown in Figure^. As expected, the number falls rapidly with 
increasing distance from the Galactic equator. In spite of the 
difference in absorption between northern and southern fields, 
there is no obvious asymmetry between counts of northern and 
southern clump giants, as expected. 



For each color, there are 4 fitted parameters: the baseline flux 
/base, the time of maximum apparent amplification to, the im- 
pact parameter at maximum amplification normalized to the 
Einstein ring radius uo = u(to), and finally, the microlensing 
event duration, i.e. the Einstein ring radius crossing time ?e = 
re/Vt for a transverse relative velocity v t . The last two quanti- 
ties depend on the Einstein radius, r 2 = 4GMD\(D s -D\)/c 2 D s , 
where M is the mass of the lens, and D\ and D s are the distances 
to the lens and source, respectively. 

While all transient flux variations can be fitted more or 
less successfully with the above expression, ideal microlens- 
ing events are characterized by achromaticity, i.e. the same 
uq, tE and to in all pass-bands. Because of the rarity of the 
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phenomenon, one generally also expects the events to be not 
repeated for a given star. A final important characteristic is 
that the detection efficiency corrected distribution of uo is flat. 
These characteristics guide the final selection of microlensing 
candidates. This selection is based on seven criteria, C1-C7, 
applied to the parameters of the microlensing fit performed on 
light curves passing the filter. The cuts we use lead to a detec- 
tion efficiency that is relatively independent of uq so the flat dis- 
tribution of mo is only moderately distorted. On the other hand, 
the cuts necessarily result in an efficiency that depends on £e 
since very short events have a significant chance of falling be- 
tween observations. The cuts are also sufficiently loose to main- 
tain a good efficiency for finding non-conventional microlens- 
ing events, e.g. those due to binary lenses. 

The first five criteria, C1-C5, are applied independently to 
the parameters of the fits for both bands and events are required 
to pass the criteria in both. CI requires that the time of maxi- 
mum apparent amplification be within the observing season so 
that the microlensing fit will yield reliable lensing parameters: 



100 



CI : 



to e observing season 



The beginning and end of the seven observing seasons were 
defined separately for each of the 66 fields by the first and last 
successful images for that field. For each field, only seasons 
with a mean sampling interval of less than 8.5 days are used. 
Sub-fields covered by CCD 2 were ignored for the first 4 sea- 
sons when the red CCD-2 was mostly out of order. 

The second cut eliminates very long candidate events: 



C2 : 



t E < 400d. 



A visual scan of events eliminated by this cut indicates that 
they are generally low amplitude events (uq > 0.8). No high 
amplitude microlensing events (uq < 0.5) were seen. Many of 
the ?e > 400 d events appear to be simple baseline shifts, with 
most occurring when the telescope optics were realigned at t — 
JD - 2,450,000 = 983. Most stars with such shifts are near 
very bright stars and it is probable that the shifts are due to 
slight changes in the wings of the PSF. 

The third cut eliminates events due to known equipment 
problems. Most importantly, we eliminate long events occur- 
ring before t = JD-2, 450, 000 = 900 by requiring 



C3 : 



t e < lOOdif f () < 900 



A visual scan of the ~ 5000 events eliminated with this cut 
indicates that almost all are due to the aforementioned baseline 
shifts at t = 983 which yielded events with fo ~ 800. We also 
eliminated 48 events on CCDs 6 and 7 between days 2323 and 
2343 that were due to an electronics problem. 

The next cut requires sufficient sampling during the event 



C4 : 



> 4 points with \t - t \ < t E 



and > 8 points with \t — to\ < 2?e . 

This cut has the largest influence on the fg dependence of the 
efficiency. 
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Fig. 6. The evolution of the uq distribution with the applied 
candidate selection criteria. The first histogram is for the 577 
events with uq < 1.0 satisfying criteria CI, C2,C3 and C4. The 
second is for the 234 events that satisfy the five criteria Cl- 
C5 and the third for the 165 events satisfying the seven criteria 
C1-C7. In the final histogram, 18 events have been eliminated 
whose light curves show additional variations outside the main 
variation or are clearly baseline shifts. The dots give the Monte- 
Carlo distribution for microlensing events satisfying all criteria. 

Cut 5 is a very loose cut on the x 1 of the microlensing fit 
calculated using only the points far from the time of maximum 
amplification, \t - fo| > 2t E : 



C5 



X base 



/AW < 10 



This cut eliminates most periodic and chaotic variable stars but 
retains single excursion events, even those that fit badly the 
simple light curve e.g. events due to binary lenses. 

Cuts 6 and 7 require reasonable agreement between the blue 
and red light curves: 



C6 



C7 



fo(blue) - ?o(red) 

f E (blue) - f E (red) 

He) 



<0.4. 



< 0.3 



For the most part, these two cuts eliminate short time-scale ex- 
cursions due to photometric problems on single exposures. 

To understand the effect of each cut it is interesting to fol- 
low the evolution of the uq distribution as successive cuts are 
applied. The first distribution in Figure|6]is for the 577 events 
with u < 1.0 that pass cuts CI, C2, C3 and C4. The distri- 
bution is dominated by events with uq > 0.5 corresponding 
to low amplitude variations. This is indicative of contamina- 
tion by low amplitude variable stars. However, already with 
just these cuts, the distribution is flat for uq < 0.4 indicating 
that there is little background for high apparent amplification 
microlensing events. 

The second distribution is for the 234 events that also sat- 
isfy C5. This eliminates most of the variable stars, and the 
mo distribution is now quite flat for mo < 0.8. Application of 
cuts C6 and C7 eliminate a few additional events, mostly near 
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event type 



number 2( ? E/f)/total 



simple 104 0.77 

strong blends (B) 5 0.02 

strong parallax (P) 4 0. 1 1 

binary non caustic (X) 2 0.02 

caustic (C) 5 0.06 

single excursion (?) 3 0.01 
chaotic/periodic variable 

or baseline shift 6 

Table 1. The tentative classifications of the 129 events satis- 
fying cuts C1-C7 and uq < 0.75. The symbols in parentheses 
(B,P,X,C,?) show how these classifications are designated in 
Figure[8]and Table[3] The last column gives the relative weight 
of each classification if they were used in the calculation of the 
optical depth using the ?e from the simple-lens fit. (In reality, 
only 120 events are used, eliminating the 3 single excursion (?) 
events and 6 chaotic/periodic variables or baseline shifts.) 



mo ~ 1, leaving 165 events with a distribution that is flat for 
m < 1.0. 

The 165 light curves were then visually examined. Most 
events fit well the simple microlensing light curve though a 
significant fraction show interesting effects. A summary of the 
scan conclusions is given in Table ^for the 129 events with 
mo < 0.75 where there is generally little ambiguity about the 
nature of the events. Non-simple microlensing events include 
four strong parallax events, i.e. those with light curves strongly 
affected by the circular motion of the Earth (P). There are 
five binary lens events with visible caustic crossings (C). Five 
events show strong source blending (B) and are discussed in 
more detail in Section|6] Two events (X) show achromatic de- 
viations from the simple lens curve that might be due either 
to non-uniform motion of the source in a binary system (i.e. 
a so-called xallarap event) or to a binary lens without caustic 
crossing. Six events show chaotic or periodic variations outside 
the primary variation and are therefore clearly not microlensing 
events. Finally three events marked (?) show single excursions 
that fit badly the simple lens curve but could conceivably be 
due to binary lenses. 

The events with mo > 0.75 are visually not very striking, 
and a non-negligible number were eliminated by scanning. We 
therefore choose to eliminate these events so that for the re- 
mainder of this paper, we consider only the 120 events selected 
by C1-C7 and classified as microlensing events with uq < 0.75, 
i.e. maximum amplifications A max > 1.6. The characteristics 
of the 120 events are listed in Table Figures |2"T1 1221 and 1231 
show the light curves of three candidate events, the first with 
the most significant improvement in x 2 over the constant-flux 
light curve, the second with the median improvement, and the 
third with the least significant improvement. All of the 120 light 
curves can be found on our web page, http://eros.in2p3.fr 

The mo distribution for the selected events is quite flat, as 
can be seen in Figure[6] It agrees with the expected distribution 
from the Monte-Carlo simulation. The fact that the observed 
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Fig. 7. The cumulative uq distribution for the 120 selected 
events. The smooth curve shows the efficiency-corrected flat 
mo distribution. 
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Fig. 8. The distribution of the reduced x 2 calculated using 
points within the range \t - ?o)| < 2?e. The 120 accepted events 
plus the 3 (?) events appear in the plot. The letters show the 
positions of the 19 non-simple events as described in Tabled 



and expected distributions fall slightly near mo = 1 is due to 
the effect of the event filter that requires a sufficient number 
of points significantly above the baseline. Figure shows the 
cumulative mo distribution for the 120 events and the efficiency- 
corrected curve for a flat uo distribution (with the observed dis- 
tribution of r E ). The Kolgomorov-Smirnov test gives a CL of 
83%. As discussed in Section|6] the sample may have a ~ 10% 
contamination at large m by blended events due to the mi- 
crolensing of faint background stars. Eliminating any randomly 
selected 10 events in the expected range 0.5 < mo < 0.75 yields 
a CL of 34%, still quite acceptable. 

The distribution of x 2 per degree of freedom calculated for 
points in the range \t - to\ < 2/g is shown in Figure [S] The 
simple events have a reasonable distribution with a mean value 
of 1 .24, a quite satisfactory value considering the difficulty of 
estimating photometric uncertainties. The other categories of 
events have, naturally, larger values of reduced x 2 , especially 
caustic events (C). 

Figure[9]shows the scatter plot of log(fE) vs. mo. The mo dis- 
tribution is, as expected, not strongly dependent on tn- Figure 
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Fig. 9. The distribution of (uo,t E ) for the 120 events passing 
all selection criteria with uq < 0.75. 

EI|shows the distribution of event magnitude and color relative 
to the center of the clump. The events show a distribution that 
is very similar to that of the ensemble of clump stars. 

The events in the sample are generally well described by 
the light curve (|8j, as indicated by the reduced x 1 distribu- 
tion (Figure [8}. Further indication that the events are indeed 
microlensing events comes from their achromaticity and time 
symmetry. Figure \^\ shows that the events have equal uq and 
?e in the blue and red as expected for microlensing events with 
negligible blending. Figure[21shows that the events have equal 
rise and fall times. This test for time symmetry was performed 
by fitting the light curves with the form © with fs in mod- 
ified by 



At = t E 



1 + a arctan 



t-t 



(10) 



The parameter a describes the asymmetry of the light curve 
with a = corresponding to a symmetric (microlensing) curve 
and a > (a < 0) corresponding to rise times less than (greater 
than) the fall time. As an example, supernova light curves gen- 
erally give a > 0.3. The observed distribution (Figure IT2t is 
symmetrically peaked at a = as expected for microlensing 
events. 

5. Optical depth and t E distribution 

To determine the optical depth and the distribution of ?e, we 
must first evaluate the detection efficiency as a function of time 
scale ?e by using Monte-Carlo simulated light curves. We su- 
perimpose artificial microlensing events, with randomly gener- 
ated parameters (impact parameter, date of maximum amplifi- 
cation and time scale), on a random sample of real light curves. 
These light curves are then subjected to the same microlensing 
search as the original light curves and we determine the fraction 
that are recovered by our detection algorithm. More precisely, 
the detection efficiency as a function of tE is given by the ratio 
between two numbers: 

- The number of events passing all selection criteria and with 
the reconstructed ?e within the considered tE bin and with 
the reconstructed uq < 0.75. 
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Fig. 10. The top panel shows the distribution of R aos relative 
to the center of the clump, /? c i llm p. The bottom panel shows the 
distribution of (B - R) aos relative to the center of the clump, 
(B - R) c \ U mp- In both cases, the dashed histogram shows the 
distribution of random clump stars. 



- The number of generated microlensing events with gen- 
erated ?e within the considered t% bin, with generated to 
within the observing period for the field in question and 
with the generated uq < 0.75. 

Figure[0]shows the efficiency as a function of ?e for a subsam- 
ple of individual fields. The efficiency is field-independent at 
the 10% level. For ?e > lOOd criterion C3 causes a slight dis- 
continuity in the efficiency for fields that were well observed 
during the first two seasons. 

The microlensing optical depth is defined as the probability 
that a given star, at a given time f, is magnified by at least 1.34, 
i.e. with an impact parameter u(t) < 1 . The optical depth is then 
given by 



71 2 != " f E,i/f( f E,i) 

2« (max) jj 



(11) 



where Af* is the number of monitored stars, Tj is the obser- 
vation period for star j, fej is the measured Einstein crossing 
time of the ith candidate and e(?E,i) is the detection efficiency. 
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Fig. 11. Tests for achromaticity of the events. The top panel 
shows the difference between red and blue minimum impact 
parameters, kq, for the 78 events with at least three points in 
the range f - f E /2 < t < to + f E /2 for both red and blue curves, 
ensuring a precise measurement of the maximum amplification 
in both colors. The bottom panel shows the difference between 
red and blue f E for the 48 events with at least three points in the 
range fo - 3f E /2 < t < to — f E /2 and three points in the range 
to + ts,/2 < t < to + 3? E /2 for both red and blue curves, ensuring 
a precise measurement of f E in both colors. 



For each event, the efficiency is taken from the corresponding 
subfield consisting of one CCD (Figure[T3l. 

The denominator of lilt is 6.62 x 10 9 star • days. Using the 
120 events in Table [5] we find a calculated optical depth aver- 
aged over all fields which is 



t= 1.68 +0.23x 10" 



(\b\) = 3.43° 



(12) 



The uncertaint y is mostly statistical , estimated following the 
prescription of Han & Gould ( 1995b) 



<t(t) = r 



< We > 



0.21 x 10 



-6 



(13) 



To the statistical error we have added in quadrature a 5% sys- 
tematic uncertainty due to blending effects, as discussed in 
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Fig. 12. Test for time symmetry of the events. The histogram is 
of the parameter a defined by JlOl i for events with fo at least 50 
days away from both the beginning and the end of the observ- 
ing season, ensuring a good measurement of a. 
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Fig. 13. Detection efficiency for individual fields as a function 
of the event time scale f E (in days) for Mo(max) = 0.75. For the 
sake of clarity, only 12 of the 66 fields are shown. The discon- 
tinuity at f E = lOOd is due to criterion C3 for fields that were 
observed during the first two seasons. 

Section[6] The location and contribution of each of the 120 can- 
didates to the measured optical depth is shown in FigurefT4l 

Because of the rapid variation of the optical depth across 
our fields, the above mean optical depth is of limited interest. 
Table [5] and Figure ^] shows the depth for various slices in 
Galactic latitude and longitude. The latitude gradient is clearly 
seen. Fitting the five points shown in Figure[T3h. one finds 



t/10" 6 = (1.62 + 0.23)exp[-a(|/>| - 3°)] 



with 



a = (0.43 ±0.16)deg" 



(14) 



(15) 



At b — -2.7 this corresponds to an optical depth gradient 
of (0.78 + 0.27) x 10 6 deg" 1 , in agreement w ith the less pre- 
cise values of MACHO JPonowski et al.l2005t) (1.06 + 0.71) x 
lO^deg" 1 and of OGT .F.-TT~( Sii mi et alJkOOfill (0.78 ± 0.84) x 
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range 


{\b\) 






(t E ) 


o-(tE) 


(fE)cor 




t/10 


-6 


1.40 < \b\ < 7.00 


3.34 


5569216 


120 


32.95 


31.07 


28.31 


±2.84 


1.68 


±0.22 


1.40 < \b\ < 2.00 


1.75 


630884 


25 


29.17 


21.85 


24.88 


±4.37 


3.52 


± 1.00 


2.00 < \b\ < 2.50 


2.26 


829123 


22 


39.29 


28.20 


30.91 


±6.01 


2.38 


± 0.72 


2.50 < \b\ < 3.00 


2.76 


976707 


24 


22.87 


15.49 


20.89 


±3.16 


1.31 


±0.38 


3.00 < \b\ < 3.50 


3.23 


931147 


25 


37.17 


37.31 


33.26 


±7.46 


2.21 


±0.62 


3.50 < \b\ < 7.00 


4.45 


2194599 


24 


36.75 


41.51 


32.60 


± 8.47 


0.92 


± 0.27 


-3.50 < b < -1.40 


2.69 


1861880 


54 


35.68 


30.22 


31.33 


±4.11 


2.42 


±0.47 


1.40 <b< 3.50 


2.50 


1512735 


42 


27.27 


23.42 


23.08 


±3.61 


1.94 


±0.42 


-6.00 < t < -3.00 


2.40 


361293 


10 


46.02 


28.13 


39.43 


±8.89 


2.90 


± 1.30 


-3.00 < e < 0.00 


2.42 


660722 


20 


29.79 


23.56 


25.65 


±5.27 


2.32 


±0.73 


0.00 <t< 3.00 


2.22 


981742 


31 


27.38 


21.75 


23.77 


±3.91 


2.20 


±0.56 


3.00 < I < 6.00 


2.53 


316650 


8 


25.93 


11.89 


23.88 


±4.20 


1.65 


±0.83 



Table 2. The measured optical depth for various slices in Galactic latitude b and Galactic longitude £. 
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Fig. 14. The location of the 120 candidates used for the calcula- 
tion of the optical depth. The area of each circle is proportional 
to the contribution of each candidate to the optical depth. 



10 6 deg l . We see no significant difference between the north- 
ern and southern fields, as expected. 

The dependence on Galactic longitude is expected to be 
much weaker. Additionally, nonuniform placement of the ob- 
serving fields makes it easy to mask a longitudinal dependence 
by the stronger latitude dependence. Figure[T5b shows the mea- 
sured optical depth as a function of longitude for EROS-2 fields 
in the latitude range 1.4° < \b\ < 3.0°. Also shown are predic- 
tions for our fields of various models. No significant longitudi- 
nal dependence is seen. 



Figure^] shows the t E distribution of the 120 events, both 
raw and efficiency corrected. The mean ?e (Table|2j is 

<? E > = (28.3 ± 2.8) d . (16) 

We see no statistically significant latitude or longitude depen- 
dence of (?e) (Table|2). The north-south difference, 

<f E >A>o - <f E >A<o = (8.2 + 5.5) d, (17) 

is not significant. Figure [n] shows the cumulative t E distribu- 
tion of the 120 events showing their contribution to the optical 
depth; 10% of the optical depth comes from the 24 events with 
f E < 13.6 d and 10% from the 4 events with t E > 120 d. 

The 15% error in the optical depth (I14> and 10% error in 
<?e) (I16> are mostly statistical, reflecting the number of events 
and the distribution of f£. We believe that it is unlikely that any 
systematic errors are this large. 

The most obvious systematic error would come from pos- 
sible contamination of the event sample with events not due 
to microlensing. This is a serious problem for microlensing 
searches in the Magellanic Clouds. We note however that the 
Magellanic microlensing rate is at least a factor 10 lower than 
in the Galactic Bulge. Furthermore, the identified Magellanic 
background events should not be a problem here since they 
concern stars that are not in the clump region (blue bumpers, 
Be stars) or background supernovae that have light curves suf- 
ficiently different from microlensing light curves to be easily 
identified for the high signal-to-noise ratio events considered 
here. At any rate, there is no indication that our sample of 120 
events is contaminated with events not due to microlensing. 
Among the indications for lack of significant contamination are 

- The^- 2 distribution (Figure[8} indicates that the function © 
is a good description of the light curves. 

- The i<o distribution is nearly flat as expected (Figures|6]and 

- The events show no sign of chromaticity (FigurelTTl. 
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Fig. 15. The top panel shows the EROS-2 measured optical 
depth as a function of Galactic latitude. The bottom panel 
shows the optical depth as a function of Galactic longitude for 
the latitude range 1.4° < \b\ < 3.0°. The filled circles are from 
this work while the o pen circle in the upp er panel is from the 
first EROS-2 analysis (lAfonso et a l. 2003b). In the upper panel, 
the solid line shows the fit dl4> . In both panels, the dotted lines 
show the prediction of the model of Bissantz & Gerhard ( 2002) 
and the dashed lines . A, B and C, show the p redictions of 
three models used bv Evans & Belokurov (2002) as described 
in Section[7J The open square in the upper panel is the predic- 
tion for the Baade Window bv lHan & Gouidlll2003l) . 



- The light curves are symmetric about the time of maximum 
amplification (Figure[QJ. 

- The optical depth has the expected dependence on the 
Galactic latitude (Figure I15> . A background of variable 
stars might be expected to yield a measured optical depth 
that is latitude-independent. 

A second systematic error could come from the contamina- 
tion of the source sample by stars not in the Galactic Bulge. The 
number of source stars that are not in fact in the Galactic Bulge 
can be estimated from the general luminosity function and stan- 
dard Galactic models. Generally, one finds that less than 1% 
of the stars in the clump box (Figure [3} could be foreground 
main sequence stars. For negative latitudes, a contribution at 
the percent level can also be expected from bright giants in 
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Fig. 16. The distribution of Einstein radius crossing times, tE- 
The top panel shows the raw distribution for the 120 events 
with mo < 0.75 with the dashed line corresponding to the 115 
events that show no strong blending. The bottom panel shows 
the distribution corrected for the ?e dependence of the detection 
efficiency. 
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Fig. 17. The normalized cumulative distribution of the contri- 
bution of the 120 candidates to the observed optical depth as 
a function of tE- 10% of the optical depth comes from the 
24 events with tE < 13.6 d and 10% from the 4 events with 
t E > 120 d. 
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the Sagittarius dwarf galaxy where stars are ~ 2.3 mag dimmer 
than low-latitude stars in the Galactic Bulge (Alard 1996). 

A third effect comes from the preferential selection of 
source stars on the near side of the Bulge because the mag- 
nitude cut on source stars (Figure [3} favors stars on the near 
side. Compared to a star at the Galactic center, stars on the 
near (far) side of the bulge are shifted upward (downward) in 
the color-magnitude diagram. Faint stars on the near side are 
moved into the accepted magnitude range while bright stars 
move out. Since there are more dim stars than bright stars, there 
is a net gain of stars on the near side. The opposite effect oc- 
curs on the far side where there is a net loss of stars. We used 
a Monte-Carlo calculation to estimate the mean position of our 
source stars and found that it is shifted towards us from the 
Galactic center by about 5% of the r.m.s. Bulge light of sight 
thickness. Since stars on the near side have a lower than aver- 
age optical depth for bulge-bulge lensing, we can expect that 
the optical depth is underestimated by of order 5%. This factor 
is not negligible compared to the statistical errors and any pre- 
cise comparison with Galactic models will require that models 
properly weight source stars by their position in the Bulge. 

A fourth effect could be the misclassification of non- 
standard microlensing events (e.g. caustic events) as variable 
stars. Our detection criteria are sufficiently liberal that several 
clear binary events were found, making up about 8% of the op- 
tical depth (Table [0. The ambiguous single excursion events 
would be only about 1% of the optical depth if they were in- 
cluded in the calculation. It therefore seems unlikely that there 
is a loss of events at the 10% level. 

A fifth systematic error could come from the use of the tz 
fitted with for events due to non-standard lenses. The prob- 
lem is minor for parallax events for which the modifications of 
are relatively small and the use of the fitted ?e makes little 
difference. On the other hand, caustic events have light curves 
that resemble (J3J only in the wings and the fitted ?e can be far 
from its true value, though for the caustic events in this study, 
the differences are only of order 10%. At any rate, the number 
of events of this category is small, so we can expect that the 
modification of the calculated optical depth should be much 
less than 10%. Thi s was the conclusion of the detailed study of 
lGiicensteinN2003l) . 

The sixth systematic effect is that due to blending, the sub- 
ject of the next section. 

6. Effect of blending on the measured optical 
depth 

Stellar blending complicates the interpretation of microlensing 
events in crowded fields. It is helpful to divide blending into 
two effects: 

- Effect I: The photometry of bright (clump giant) stars is af- 
fected by the background of faint stars that are randomly 
placed with respect to the bright star. In the photometry, the 
sky background is assumed to be a smooth function so if 
the number of faint stars inside the seeing disk of the bright 
star is greater than the average number, the baseline flux of 
the bright star will be overestimated and the amplification 



due to microlensing underestimated. The opposite will oc- 
cur if the number of faint stars is smaller than average. We 
can therefore expect that on average amplifications are cor- 
rectly estimated. However, as emphasized in the previous 
section, the magnitude cut (Figure used to select source 
stars favors stars whose baseline fluxes are overestimated. 
We therefore expect a slight underestimation of the ampli- 
fication. 

- Effect II: The microlensing of a faint star within the seeing 
disk of a bright star will cause the reconstructed flux of the 
bright star to vary in time, yielding an apparent microlens- 
ing of the bright star. 

By itself, effect II will clearly cause one to overestimate 
the optical depth if (II It is used since some of the events will 
not be due to the clump giants. On the other hand, we ex- 
pect that effect I goes in the opposite direction since amplifica- 
tions of clump giants are on average underestimated. Because 
of this, the reconstructed ?e found using (|8} are systematically 
underestimated while the reconstructed uq are overestimated. 
In the formula ( II It for the optical depth, the terms in the sum 
are therefore underestimated on average. Furthermore, fewer 
events are included in the sum because the overestimation of uq 
causes some events to migrate beyond uo{max). Effect I there- 
fore, by itself, causes the optical depth to be underestimated. 

To estimate the importance of blending we have performed 
two types of tests with artificial stars placed on CCD images. 
The first type uses totally synthetic images in which all stars 
are artificial. This was used to quantify the sum of effects I and 
II. The second type uses real EROS-2 images in which artificial 
clump giants are added at random places on the images. This 
was used to quantify effect I. 

The production and analysis of the totally synthetic im- 
ages were described i n the first EROS-2 Galactic Bulge article 
dAfonso et al] Eo03b). Series of CCD images were produced 
by placing on them stars taken randomly from the luminosity 
distributi on observed for th e Baade window in HST deep im- 
ages (Holtzm an et al.lll9 98). The HST luminosity distribution 
extends to 9 magnitudes fainter than clump giants and there- 
fore allowed us to simulate the bright stars used in the EROS-2 
analysis as well as the fluctuating background of faint stars. 
Microlensing events were assigned to a subsample of stars and 
the resulting flux change was taken into account when placing 
charge on the CCD. The sequence of CCD images were then 
reconstructed in the same way as normal EROS-2 images and 
events were searched for. 

The results of this study indicated that reconstructed ampli- 
fications are indeed underestimated for clump giant stars lead- 
ing to an underestimation of the optical depth (effect I) but that 
this is compensated to a precision of about 5% by events due 
to backg round stars (effect II). For the analysis of lAfonso et alJ 
(2003b), about 18% of the recovered optical depth was due to 
these faint stars compensating for the 15% underestimation for 
those events due to bright stars. For the analysis presented here, 
the effects are smaller because brighter stars are used and the 
amplification threshold is higher (1.6 instead of 1.34). We now 
find that only about 6% of the recovered optical depth is ex- 
pected to be due to faint background stars. 
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Fig. 18. The ratio of recovered and generated values of f E for 
artificial clump giants placed on EROS-2 field 610. 



Our second test used partially synthetic images on which 
we placed artificial clump giants on the real EROS-2 CCD im- 
ages. As in the first test, microlensing events were simulated 
by giving the artificial giants time-dependent fluxes. These so- 
called spiked images were then photometered by the normal 
procedure and the events recovered and fitted with the mi- 
crolensing light curve. 

Figure ED shows the comparison between recovered and 
simulated f E . This parameter is correctly estimated on average 
with 85% of events yielding a f E within 20% of the generated 
f E . The 15% of events in the wings are due to stars whose base- 
line flux was strongly overestimated or underestimated because 
of blending. The nearly symmetric shape of the distribution in 
Figure^] means that a bias in the reconstruction of f E results 
mostly from the fact that the magnitude cuts favor stars whose 
flux was overestimated. The distribution in Figure ^] and the 
associated correlation with reconstructed baseline flux was thus 
combined with the observed magnitude distribution to make a 
Monte-Carlo estimate of the optical depth bias. The results con- 
firm that effect II causes ~ 5% underestimation of the optical 
depth. 

Since these images only give information on effect II, we 
supplement it with a simple numerical calculation to estimate 
the contribution of faint stars in the seeing disk. A microlensed 
faint star superimposed on a bright star yields a light curve 
given by 



F(t) = F s 



(1-/) +/ 



u 2 + 2 



u iu 2 +4 



(18) 



where / is the ratio between the faint star's baseline flux and 
the total reconstructed baseline flux. For / < 1, the recon- 
structed amplification is much less than the real amplification 
yielding reconstructed uq and f E respectively much greater than 
and much less than the real uq and f E . In practice, this means 
that stars more than 4 magnitudes fainter than clump giants 
yield no observable events because their reconstructed f E are 
less than 5 d. The HST luminosity function gives on average 
about 3 stars within 2 arcsec of a clump giant and within 4 mag- 
nitudes of the clump magnitude. Integrating the HST function, 



we can expect that about 10% of clump giant events will be 
due to faint background stars. Since the events are much shorter 
than real clump giant events, they contribute only about 4% to 
the measured optical depth, not far from the above estimates. 

Because of the approximate compensation, we make no 
blending correction to the optical depth. Instead, because of 
the difficulty in estimating the effect precisely, we choose to in- 
clude a 5% uncertainty in the optical depth, i.e. an uncertainty 
equal to the estimated size of each of the compensating effects. 
This uncertainty is, at any rate, less than the statistical uncer- 
tainty of the present measurements. Future high statistics mea- 
surements of the optical depth will certainly have to be much 
more careful about blending, especially if faint source stars are 
used. 

Some strongly blended events can be found by fitting 
light curves to a microlensing curve Jl 81 . Since all events are 
blended to some extent, the number of detected blends is a 
sensitive function of the photometric precision. Five EROS-2 
events showed a significant improvement for this curve over 
that with the simple curve l|8}. The events are characterized by 
excess amplification in their wings due to the fact that the real 
f E is greater than the f E fitted with 

The absence of strong blending in most events is seen in 
Figure 119b showing f E (blend) vs. fe(no - blend). While the 
cloud of points is very wide because of the degeneracies in 
the fit, for the most part there is no systematic difference be- 
tween the two f E , whereas f E (blend) > f E (no - blend) for truly 
blended events. This is confirmed in Figure [T^b. which shows 
the logarithm of the ratio of the two f E 's both for the observed 
events and for the simulated clump giants. The observed distri- 
bution is very similar to the Monte-Carlo distribution. 

The absence of strong blending for our events is not incon- 
sistent with the results of the OGLE-II group who observed 
a large number of strongly blended events on clump giants. 
Most of the OGLE-II blended events are very low amplifica- 
tion events, («o(no blend) > 1) with only 4 of their 26 events 
with Mo(no blend) < 0.75 being characterized as blends. This 
proportion of blends is somewhat higher than our estimates, a 
fact that may be due to their superior photometric precision. At 
any rate, the number is sufficiently small to be of minor impor- 
tance for optical depth determinations with uq < 0.75. 

7. Discussion and Conclusions 

The EROS-2 optical depth measurement presented here is in 
remarkably good agreement with the other clump-giant mea- 
surem ents shown in Figure l20l i.e. those of the MACHO 
group JPopowski et al.l Eo05) and OGLE-II group dSumi et alJ 
2006). O nly the ori ginal low statistics MACHO measurement 
dAlcocketalJl997h is a bit too far from our fitted optical depth 

<HS. 

Figure [21 shows that our optical depth measurements are 
in reasonable agreement with some recent calculations based 
on models of the Galactic Bar and Disk. The dotted line in 
Figure [15] is the pr e dictio n of the non-parametric model of 
Bissantz & Gerhard (2002). The strong latitude dependence 
is correctly predicted as is the nearly flat longitude depen- 
dence. The three dashed lines correspond to calculations of 
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Fig. 19. The fitted f E assuming blending (eauationll8>. vs. the 
fitted value assuming no blending (equation[S}. The filled cir- 
cles correspond to the five events showing strong blending. The 
bottom panel shows the logarithm of the ratio of the two fg for 
simulated events (solid line) and observed events (dashed line). 
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Fig. 20. Optical depth as a function of Galactic latitude. The 
filled circles are t his work while the open circle is the first 
EROS-2 analysis ( Afo nso et al.ll2003bl) . The solid line shows 
the fit (1141 . T he measurements of th e MACHO group are the 
open squares (Popowski et al. 2005) and the crossed square of 
the early analysis of Alcock et al. ( 1997). The crosses are from 
the OGLE-II group isumi et all200d) . 



lEvans & Belokurovl (12002) based on the three models of the 
inner G alaxy. Curve A is based on t he model of iBinnev et alJ 
dl997l) . curve B on the mod el of of iDwek et all i ll 9951) . and 
curve C on the model of of Freudenreich ( 1998). (For these 
three models, we show the predictions without corrections due 
to spiral structure; these would give at most a 20% increase to 
the calculated optical depth.) All three models have a common 
total mass of 1 .5x 1O 1O M within 2.5 kpc of the Galactic Center. 
The data in Fi gure LL5Jindicate that the "swollen bar" model of 
iFreudenreichl dl998l) is not favored by the data, unless the as- 
sumed mass is reduced by ~ 30%. 

Finally, Figur e \15\ shows the Baade window prediction of 
lHan & G ould (2003). The calculated optical depth is within 
one standard deviation of the curve given by (I14> . 

We make no attempt here to separate the Disk and 
Bulge Q contributions so ou r concl usion is simply that the 

model of Bissantz & Gerhard (2002) or models A and B of 

I ^ — 

Evans & Belokurov ( 2002) give a good estimate of the total op- 
tical depth ([0. We note that these models place the bulk of the 
mass in normal stars and gas so the agreement with our mea- 
surements indicate that there is no need for additional mass in 
a non-lensing form. In models with cold dark matter, the dark 
component may make a significant contribution to the mass in 
the inner parts of the Milky Way. F or example, in the models of 
iKlvpin. Zhao & Somerville1(l2002l) . the cold dark matter contri- 
bution to the rotation curve at 1 kpc from the Galactic Center 
is between 30% and 50%. Such a contribution would, by it- 
self, lower the optical depth by roughly the same factor. Our 
optical depth measurement, with its 15% uncertainty, therefore 
provides a useful constraint on such models. 

Galactic models generally suppose a smooth distribution 
of lenses and therefore predict a smooth longitude and lat- 
itude dependen ce of the optical depth. The MACHO group 
(Popows ki et alJl2005l) observed 9 events near ({ = 2.9°, b = 
-2.9°) (their field 104), yielding an optical depth about 2 stan- 
dard deviations above the expectation based on their other 
fields. As seen in Figure[2] we do not see an excess in this di- 
rection. Of their 9 events, 4 events occurred after the beginning 
of EROS-2. Of these 4 events, two are not in an Eros field and 
two are in our field 611. Both of these events were found (in 
Tabled MACHO events 104.205 15.498 and 104.20640.8423). 
No other EROS events were found in this region in the three 
bulge seasons after the shutdown of MACHO. The fact that 
EROS-2 saw only two events for a solid angle and time period 
equivalent to that in which MACHO found 9 events suggests 
that the MACHO excess is a statistical fluctuation. 

While our optical depth measurement is in good agreement 
with other measurements and with Galactic models, agreement 
on the tE distribution is less satisfactory. Our value of (f E ) is 
in excelle nt agreement with the value, (28.1 + 4.3) d, found by 
OGLE-II JSumi et alJ2006h . The EROS-2 and OGLE-II values 
are, however, significantly higher than that calculated using the 
62 MACHO events (Table 3 of dPopowski et al]l2005h : <f E > = 
(2 1 .6 ±3) d. The discrepancy lies entirely in the large number of 
MACHO events with f E < 10 d. The number of MACHO events 
with (r E < 5, 5 < t E < 10, t E > lOd) is (6, 14,42) whereas the 
corresponding numbers for EROS-2 are (0, 10, 1 10). Part of the 
differences in event numbers is due to the fact the the MACHO 
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efficiency is relatively greater than the EROS-2 efficiency for 
t E < 10 d. Using the 14 MACHO events with 5d < t E < 10 d 
and the relative efficiencies of the two experiments, we estimate 
that EROS-2 should see 24 ± 7 events in this range compared 
to the 10 that are seen. The very low EROS-2 efficiency for 
?e < 5 d makes it unsurprising that we found no events in this 
range. 

Whatever the source of the discrepancy, these short events 
have relatively little effect on the measured optical depths. 
The events with t E < 10 d constitute only about 12% of the 
MACHO optical depth and only 4% of the EROS-2 optical 
depth. 

The Galactic models are capable of reproducing the dis- 
tribution by adopting an appropri ate stellar mass function. The 
model of IWood and Maol (12005) is in good agreement with 
our distribution an d with that of OGL E-II, both with (t E ) ~ 
28 d. The model of lBissantz et alJ J2004I) is in good agreement 
with the dist ribution reported by M ACHO with <f E ) ~ 21 d. 
According to lWood and Maol J2005h . the different (t E ) predic- 
tions of the two models is mostly due to different adopted 
mass functions l eading to a mean lensi ng mass of 0.35Mo fo r 
IWood and Maol ( 120051) and 0.1 1M for lBissantz et alJ j2004). 
Kinematical differences between the two models may also have 
a non-negligible effect. 

The present comparison with models is still primarily lim- 
ited by the low number of observed events. Since the EROS-2 
survey uses most of the fields in the Bulge region that are eas- 
ily observed in optical bands, a large increase in the number of 
monitored clump giants would require infrared observati ons o f 
highly obscured regions near the Galactic Center (Gould 1995). 
Of course the number of events can also be increased by using 
dim stars but, in this case, reliable results would require a good 
understanding of blending. 
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event 


r.a 


dec. 


b 


I 


tE 




to 




note 



607-1-4-11480 
607-3-3-9841 
607-5-1- 14742 
607-5-1-29607 
607-5-4- 20285 
607-6-3- 29248 
610-0-1- 17627 

610- 2-3-4703 

611- 1-4-24874 

611- 3-4- 18476 

612- 7-4- 8538 

613- 5-1-23714 
615-1-2- 22556 
619-0-4- 7664 
619-5-4- 19755 
619-7-4- 14204 
624-4-3- 3703 
627-7-2- 14994 

13467 
18520 
5329 
30438 
20157 
2992 
27217 
20099 
26926 
15890 
5389 
2508 
■9143 
8464 
■7864 
■20818 
19218 
8951 



2-0-1- 
2-0-4- 
2-1-1- 
2-3-3- 

2- 5-1- 

3- 4-3- 
3-4-3- 
3-5-2- 
3-5-3- 

3- 6-4- 

4- 3-2- 
4-3-3- 
4-4-4- 
4-5-2- 
4-7-1- 
4-7-1- 

4- 7-2- 

5- 4-1- 



5-4-3-22912 
5-4-4- 11338 



5- 5-2- 

6- 1-3- 
6-5-3- 
8-3-1- 



19037 
■21678 
18380 
1345 
8-4-2-21125 
8-4-2-21623 

8- 7-3- 16227 

9- 7-4-7143 

30- 2-3- 19447 

31- 3-3-35856 
61-0-3- 8256 
70-3-1- 18517 
70-3-1- 18797 
70-3-4- 28566 
70-4-1- 17316 
70-5-2- 24247 
70-6-2-31993 



270.2603 
270.2787 
270.0417 
270.0288 
270.2924 
269.8908 
270.4365 
270.6194 
271.0400 
271.1397 
271.8072 
271.5512 
271.8863 
272.2137 
272.7316 
272.7335 
273.1066 
273.4693 
267.8171 
268.1808 
268.3954 
268.4314 
268.2981 
268.1774 
268.2445 
268.3380 
268.6611 
268.2890 
269.2294 
269.4168 
268.9565 
269.0497 
269.1321 
269.1327 
269.2468 
268.8382 
269.0256 
268.9376 
269.2080 
270.0438 
270.1206 
270.1013 
269.7872 
269.6680 
270.2820 
271.0161 
273.9211 
274.4120 
274.4475 
265.4297 
265.3838 
265.5989 
265.0553 
265.3672 
264.9205 



-29.0032 
-29.3148 
-29.5630 
-29.6350 
-29.8155 
-30.0002 
-28.9678 
-29.2331 
-27.7431 
-28.1256 
-32.0593 
-30.1737 
-26.7402 
-28.4529 
-29.0649 
-29.4067 
-30.1842 
-26.4997 
-30.2087 
-30.4511 
-30.1464 
-30.6361 
-30.9266 
-29.4347 
-29.5728 
-29.6770 
-29.5536 
-30.0199 
-30.2536 
-30.0798 
-30.6300 
-30.6251 
-30.8306 
-30.9075 
-31.0564 
-29.0702 
-29.1386 
-29.2330 
-29.2700 
-30.3994 
-31.0448 
-27.7379 
-28.4151 
-28.3593 
-28.5282 
-31.5437 
-28.5940 
-27.2411 
-22.4102 
-34.2263 
-34.2285 
-34.4414 
-34.5743 
-34.7616 
-35.1415 



-2.942 
-3.110 
-3.053 
-3.079 
-3.367 
-3.156 
-3.059 
-3.328 
-2.920 
-3.184 
-5.593 
-4.493 
-3.090 
-4.170 
-4.862 
-5.026 
-5.679 
-4.217 
-1.710 
-2.103 
-2.109 
-2.383 
-2.432 
-1.585 
-1.705 
-1.828 
-2.008 
-1.965 
-2.786 
-2.839 
-2.770 
-2.837 
-3.002 
-3.041 
-3.200 
-1.897 
-2.073 
-2.055 
-2.277 
-3.468 
-3.844 
-2.195 
-2.290 
-2.171 
-2.724 
-4.755 
-5.557 
-5.308 
-3.067 
-2.068 
-2.036 
-2.300 
-1.989 
-2.306 
-2.197 



1.698 

1.441 

1.121 

1.054 

1.020 

0.680 

1.809 

1.666 

3.132 

2.849 

-0.195 

1.288 

4.378 

3.057 

2.769 

2.478 

1.984 

5.306 

-0.441 

-0.482 

-0.125 

-0.524 

-0.831 

0.381 

0.294 

0.248 

0.501 

-0.066 

0.164 

0.397 

-0.280 

-0.233 

-0.370 

-0.435 

-0.509 

0.993 

1.020 

0.900 

0.992 

0.411 

-0.102 

2.709 

1.986 

1.979 

2.114 

-0.118 

3.712 

5.095 

9.292 

-4.914 

-4.936 

-5.016 

-5.374 

-5.388 

-5.908 



12.780 

17.354 

18.181 

61.540 

36.233 

15.967 

18.775 

10.458 

27.429 

13.708 

70.213 

20.356 

14.728 

9.705 

16.290 

116.099 

17.672 

43.505 

59.171 

22.388 

26.760 

38.098 

97.853 

65.395 

14.034 

14.104 

23.762 

38.201 

16.954 

13.941 

21.513 

17.209 

137.482 

14.483 

37.820 

10.424 

30.119 

7.704 

7.912 

13.550 

23.286 

78.728 

34.932 

43.148 

19.994 

50.982 

23.998 

13.824 

37.373 

90.589 

91.418 

29.865 

7.140 

71.090 

46.872 



0.689 
0.691 
0.067 
0.133 
0.498 
0.709 
0.347 
0.420 
0.619 
0.699 
0.233 
0.175 
0.126 
0.408 
0.507 
0.264 
0.635 
0.041 
0.519 
0.239 
0.506 
0.552 
0.088 
0.476 
0.129 
0.645 
0.437 
0.066 
0.527 
0.233 
0.044 
0.456 
0.728 
0.614 
0.688 
0.724 
0.492 
0.640 
0.000 
0.060 
0.208 
0.726 
0.358 
0.124 
0.469 
0.510 
0.398 
0.318 
0.654 
0.181 
0.207 
0.425 
0.567 
0.061 
0.246 



1274.70 

716.60 

2113.40 

636.60 

1315.50 

2137.00 

1775.80 

1988.40 

683.50 

998.40 

1019.10 

2144.10 

1599.30 

523.40 

1732.70 

1254.70 

1414.90 

1958.20 

2383.30 

1707.10 

2122.60 

996.30 

1726.80 

2100.30 

1685.60 

1759.70 

1257.30 

955.90 

1283.20 

1411.40 

1719.10 

1985.60 

2086.20 

2180.60 

1833.10 

1093.70 

2063.00 

1271.10 

1651.60 

1763.90 

2077.50 

1714.00 

948.70 

2129.90 

1301.90 

1450.80 

1698.10 

1628.90 

1743.70 

1707.40 

2400.50 

1602.40 

1428.90 

1842.60 

1456.10 



0.576 
0.577 
0.713 
0.823 
0.810 
0.560 
0.671 
0.497 
0.798 
0.706 
0.867 
0.722 
0.650 
0.474 
0.608 
0.677 
0.588 
0.692 
0.793 
0.682 
0.783 
0.908 
0.863 
0.782 
0.582 
0.554 
0.706 
0.489 
0.822 
0.753 
0.797 
0.744 
0.930 
0.707 
0.912 
0.537 
0.825 
0.411 
0.411 
0.625 
0.688 
0.907 
0.801 
0.887 
0.803 
0.914 
0.758 
0.503 
0.781 
0.964 
0.961 
0.871 
0.330 
0.890 
0.587 



SC30-57488 
M-l 13.19192.365 

M-118. 18797.1397 
M-118. 19184.939 

SC30-636963 

M-l 04.205 15. 498 
M-l 04.20640. 8423 



M-l 10.22455.842 
P 



SC37-645044 
SC3-371229 
SC4-522952 
SC3-91382 



B 
P 

SC22-414328 

B 

C 

XM-40 1.48408 .649 
SC30-165305 
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71-3-1-21364 266.2818 -33.3554 -2.214 -3.797 47.100 0.073 892.10 0.911 

71- 7-1-28518 266.4633 -34.0936 -2.727 -4.332 36.263 0.177 2024.70 0.906 C 

72- 1-1-22000 266.3057 -34.4337 -2.792 -4.688 18.317 0.131 2343.10 0.722 

72- 4-2- 12551 265.8097 -35.1914 -2.841 -5.544 21.568 0.461 1442.30 0.804 

73- 3-4-7001 267.3455 -33.4813 -3.040 -3.418 14.286 0.043 896.30 0.739 

74- 5-3- 16791 267.4200 -35.0934 -3.917 -4.730 23.496 0.180 2177.00 0.793 

76- 1-3- 12658 268.3605 -31.6522 -2.845 -1.417 82.468 0.676 1046.30 0.900 

77- 0-1-20761 267.5647 -33.1794 -3.042 -3.065 18.932 0.487 928.60 0.711 
77-0-3-7242 267.7949 -33.0623 -3.149 -2.862 128.844 0.151 1628.20 0.883 P 
77-1-4- 10694 268.3326 -33.2374 -3.625 -2.765 24.832 0.419 2005.60 0.853 

77- 4-1- 13517 267.7165 -33.8108 -3.473 -3.525 23.493 0.144 1426.00 0.856 

78- 4-1-3363 268.7127 -31.8944 -3.225 -1.461 8.176 0.430 1387.10 0.492 

79- 2-1-8900 268.4667 -33.0348 -3.620 -2.534 29.763 0.344 1680.00 0.832 

80- 3-2-7178 269.8741 -32.2297 -4.246 -1.216 35.887 0.049 1047.00 0.876 

80- 4-2- 10880 269.4209 -32.6063 -4.100 -1.740 10.348 0.733 2040.10 0.548 

81- 3-2-24224 269.7945 -33.7675 -4.946 -2.545 30.017 0.253 1341.00 0.887 

82- 2-3-24164 270.3955 -32.2755 -4.652 -1.018 13.832 0.285 2357.80 0.631 

83- 1-1-8206 270.5674 -33.2471 -5.253 -1.758 13.159 0.180 1315.10 0.614 

84- 4-1- 15036 271.1088 -33.0547 -5.556 -1.351 19.032 0.000 2176.00 0.727 

101- 1-1-20136 263.7159 -29.3984 1.743 -1.622 10.623 0.094 2351.60 0.292 B 

102- 0-3- 17623 262.5237 -29.4664 2.577 -2.223 25.501 0.152 1322.90 0.593 
102-3-1- 13451 262.6484 -29.7646 2.322 -2.420 12.616 0.290 1093.50 0.554 
102-5-2- 15625 262.7204 -30.2758 1.990 -2.822 24.927 0.139 2197.60 0.587 

102- 5-3-29920 262.9454 -30.2184 1.858 -2.671 33.810 0.586 1348.90 0.642 

103- 1-4-21106 263.1299 -28.1925 2.829 -0.858 23.528 0.310 2380.90 0.630 

103- 3-4-27325 263.1901 -28.5500 2.590 -1.137 10.030 0.662 1777.20 0.463 

104- 3-3-29860 262.2052 -29.1929 2.959 -2.134 12.518 0.142 2064.20 0.572 

104- 3-4-27383 262.1975 -29.3197 2.895 -2.246 20.467 0.575 1984.30 0.651 

105- 5-1- 14483 261.1343 -29.5197 3.553 -2.900 7.628 0.682 1636.30 0.334 

105- 7-4- 12278 261.2801 -30.0158 3.171 -3.256 25.769 0.356 1604.50 0.739 

106- 1-4- 11177 265.1363 -27.0417 1.947 1.043 30.726 0.324 1624.90 0.565 
106-4-1-8809 264.7175 -27.5637 1.985 0.404 5.775 0.298 1701.50 0.291 
106-4-4-24841 264.8552 -27.8152 1.748 0.251 24.738 0.453 1073.50 0.785 B 
106-4-4-31028 264.8404 -27.8541 1.739 0.211 17.878 0.535 1623.80 0.687 
106-5-1- 15482 265.0711 -27.6027 1.699 0.531 19.899 0.122 1293.10 0.610 
106-5-1-26077 264.9651 -27.6713 1.742 0.424 11.974 0.593 2099.10 0.462 
106-5-4-24760 265.1471 -27.8188 1.527 0.381 25.477 0.093 2524.00 0.679 

106- 6-1- 13061 264.6275 -27.9500 1.847 0.032 36.059 0.283 1687.00 0.508 C 

107- 0-2-21213 263.4739 -26.7017 3.382 0.578 124.365 0.025 1250.10 0.776 
107-1-2-24707 263.7571 -26.7073 3.166 0.702 16.090 0.404 2080.40 0.543 
107-3-3-31150 263.9854 -26.9224 2.878 0.621 19.904 0.591 1309.50 0.754 
107-4-1-8762 263.3894 -27.1358 3.210 0.166 35.398 0.489 1754.90 0.758 
107-4-4-23806 263.7192 -27.3873 2.827 0.100 53.574 0.650 1055.00 0.881 
107-5-2-26345 263.8999 -27.3882 2.692 0.182 14.064 0.197 1710.70 0.576 
107-7-1-22583 263.8056 -27.5884 2.654 -0.033 8.574 0.075 1635.90 0.415 B 

107- 7-2- 167 263.9081 -27.6208 2.561 -0.014 9.466 0.220 1616.40 0.455 

108- 3-1- 12859 266.2558 -25.8209 1.736 2.604 12.864 0.499 1658.60 0.585 
108-3-1-27549 266.2520 -25.9237 1.685 2.513 11.375 0.197 888.20 0.504 
108-4-1- 13544 265.9802 -26.1798 1.760 2.169 53.105 0.689 1322.70 0.854 

108- 4-1-23818 265.7825 -26.2489 1.874 2.019 99.504 0.672 2532.60 0.823 

109- 4-3- 15519 265.2998 -25.5322 2.621 2.417 17.836 0.417 1365.00 0.686 
112-0-2- 12627 266.7715 -23.3296 2.629 4.995 21.280 0.615 1317.10 0.639 
112-2-2-25953 266.8246 -23.7485 2.371 4.655 20.012 0.700 1648.60 0.765 
112-5-2-24472 267.1992 -24.0774 1.907 4.541 23.119 0.676 1011.60 0.558 
112-6-1-4305 266.7831 -24.1610 2.190 4.279 10.871 0.293 2473.80 0.363 

112- 7-3- 14923 267.3655 -24.2420 1.693 4.474 54.797 0.612 1813.50 0.789 

113- 1-2-28239 266.3130 -23.4511 2.926 4.681 23.803 0.066 1967.50 0.768 
113-7-1-23134 266.3956 -24.3175 2.411 3.967 26.167 0.294 863.60 0.817 
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114-3-1- 


11753 


265.5659 


-23.5507 


3.458 


4.254 


18.387 


0.076 


2044.50 


0.728 




114-5-1- 


16265 


265.5462 


-23.9263 


3.276 


3.919 


67.940 


0.548 


1664.30 


0.810 


X 


117-6-3- 


13250 


268.0484 


-22.5049 


2.042 


6.293 


6.270 


0.082 


1643.30 


0.203 




117-7-1- 


12630 


268.2245 


-22.5131 


1.898 


6.366 


24.160 


0.100 


1794.70 


0.761 




124-1-1- 


16687 


259.9039 


-27.8925 


5.362 


-2.066 


19.173 


0.180 


1410.20 


0.688 




124-7-2- 


20832 


259.9118 


-29.1101 


4.664 


-3.106 


201.882 


0.311 


2075.30 


0.854 


P 


124-7-3- 


21344 


260.0548 


-28.9748 


4.638 


-2.925 


47.000 


0.295 


918.50 


0.851 





Table 3: The 120 events used for the measurement of the optical depth. 
The first column gives the identification: field, CCD, quadrant, star num- 
ber. The time of maximum amplification, to, is shown in the eighth col- 
umn as (JD-2,450,000). The efficiency, e, used in the optical depth cal- 
culation is shown in the ninth column. The final column signals events of 
special character: P (parallax), X (xallarap or binary lens without caus- 
tic), C (caustic binary lens), B (strong blend). Eve nt numbers in the final 
column preceded by M refer to MACHO events dPopowski et alJl2 005 ) 
and those preceded by SC refer to OGLE-II events ISumi et al.l (12006). 
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Fig. 21. The light curves of EROS-2 microlensing candidate 6-5-3-18380. Of the 120 candidates, this one has the largest value 
of^ 2 (flat) - x 1 (microlensing) = 221600. The two top panels show R ems and B e ros as a function of time (JD-2,450,000). The two 
bottom panels show a zoom for the time interval fo - 2?e < f < fo + 2fE. 
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Fig. 22. The light curves of the EROS 2 microlensing candidate 107-3-3-31150. Of the 120 candidates, this one has the median 
value of^- 2 (flat) - x 1 (microlensing) = 6810. The two top panels show R ao& and B e ros as a function of time (JD-2,450,000). The 
two bottom panels show a zoom for the time interval to - 2?e < t < to + 2?e- On the bottom two panels, the solid line is the 
blend-free simultaneous fit for both colors, the dashed line the blend-free single-color fit, and the dotted line the simultaneous fit 
with blending. 
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Fig. 23. The light curves of the EROS 2 microlensing candidate 70-4-1-14316. Of the 120 candidates, this one has the smallest 
value of^ 2 (flat) - ^(microlensing) = 399. The two top panels show R el0S and Z? e ros as a function of time (JD-2,450,000). The 
two bottom panels show a zoom for the time interval to - 2t% < f < fo + 2?e- On the bottom two panels, the solid line is the 
blend-free simultaneous fit for both colors, the dashed line the blend-free single-color fit, and the dotted line the simultaneous fit 
with blending. 



